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an extended Higgs sector, a very light CP-odd Higgs boson can naturally arise making it possible 
for a very light neutralino to annihilate efficiently enough to avoid being overproduced in the 
early Universe. 

In this article, we explore the characteristics of a supersymmetric model needed to include a 
very light neutralino, 100 MeV < m~o < 20 GeV, using the NMSSM as a prototype. We discuss 
the most important constraints from Upsilon decays, h S7, fi^fi" and the magnetic 

moment of the muon, and find that a light bino or singlino neutralino is allowed, and can be 
generated with the appropriate relic density. 

It has previously been shown that the positive detection of dark matter claimed by the 
DAMA collaboration can be reconciled with other direct dark matter experiments such as 
CDMS II if the dark matter particle is rather light, between about 6 and 9 GeV. A singlino or 
bino-like neutralino could easily fall within this range of masses within the NMSSM. Addition- 
ally, models with sub-GeV neutralinos may be capable of generating the 511 keV gamma-ray 
emission observed from the galactic bulge by the INTEGRAL/SPI experiment. 

We also point out measurements which can be performed immediately at CLEO, BaBar and 
Belle using existing data to discover or significantly constrain this scenario. 
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I. INTRODUCTION 



Despite the substantial effort which has gone into its detection, the nature of dark matter 
remains unknown [1]. The dark matter candidates which have received the most attention fall 
into the category of Weakly Interacting Massive Particles (WIMPs), which can emerge from a 
variety of theoretical frameworks, including supersymmetry. Of the supersymmetric candidates 
for dark matter, the lightest neutralino is often considered to be the most attractive. 

Neutralinos produced in the early Universe must annihilate into Standard Model particles 
at a sufficient rate to avoid overproducing the density of dark matter. Within the framework of 
the Minimal Supersymmetric Standard Model (MSSM), the lightest neutralino can annihilate 
through a variety of channels, exchanging other sparticles, Z bosons, or Higgs bosons. The 
masses of sparticles such as sleptons or squarks, as well as the masses of Higgs bosons, are limited 
by collider constraints, with typical lower limits of around ~100 GeV. For lighter neutralinos, it 
becomes increasingly difficult for these heavy propagators to generate neutralino annihilation 
cross sections that are large enough. The most efficient annihilation channel for very light 
neutralinos in the MSSM is the s-channel exchange of a pseudoscalar Higgs boson. It has been 
shown that this channel can, in principle, be sufficiently efficient to allow for neutralinos as light 
as 6 GeV [2]. Such models require a careful matching of a number of independent parameters, 
however, making viable models with neutralinos lighter than ~20 GeV rather unlikely [4]. 
Measurements of rare 5-decays are also particularly constraining in this regime. If we do not 
require that the LSP be the dominant component of dark matter, its mass can be zero [5]. 

More generally speaking, Lee and Weinberg have demonstrated that a fermionic dark matter 
candidate which annihilates through its couplings to the weak gauge bosons must be heavier 
than a few GeV to avoid over-closing the Universe [6]. Therefore, if a neutralino is to be very 
light, it requires another annihilation channel which enables it to sufficiently annihilate in the 
early Universe. This can be provided within the context of the Next-to-Minimal Supersym- 
metric Standard Model (NMSSM) by the lightest of the two CP-odd Higgs bosons, which can 
be considerably lighter than the single CP-odd Higgs boson of the MSSM without violating 
collider constraints. Furthermore, it has been shown that models within the NMSSM which 
require the smallest degree of fine tuning often contain a light CP-odd Higgs boson [7]. 

In addition to these theoretical arguments, there are experimental motivations to consider 
light dark matter particles. The observation of 511 keV gamma rays from the galactic bulge [8] 
indicates the presence of a Gaussian profile of low-velocity positrons throughout our galaxy's 
inner kiloparsec. It is challenging to explain this observation with traditional astrophysics [9]. 
Annihilating [10] or decaying [11] dark matter particles have been suggested as a possible source 
of these positrons. If such a dark matter particle were in the mass range usually considered, 
however, their annihilation would produce positrons with far too much energy to annihilate at 
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rest. Furthermore, they would almost certainly generate far too many gamma-rays and violate 
the constraints placed by EGRET [12]. Thus a dark matter candidate capable of generating 
the observed 511 keV line must be exceptionally light. 

Additionally, it has been shown that the claims of dark matter detection made by the DAMA 
collaboration [13] can be reconciled with null results of CDMS II [14] and other experiments if 
one considers a WIMP lighter than approximately 10 GeV [15, 16]. 

In this article, we explore the phenomenology of supersymmetric models with a neutralino 
in the mass range of 100 MeV to 20 GeV within the context of the NMSSM. We find that 
many such models can be found which are not highly fine tuned and are consistent with all 
constraints including direct collider searches, rare decays, and relic abundance considerations. 
We find examples of consistent models in which a light neutralino can potentially produce the 
511 keV emission observed by INTEGRAL as well as models that can potentially reconcile 
DAMA with CDMS II. However, we have not found models in which all these observations can 
be simultaneously explained. 

II. NEUTRALINO DARK MATTER WITH A SINGLET HIGGS 

The simplest possible extension of the particle content of the MSSM is the addition of 
a new gauge singlet chiral supermultiplet. There are several ways to do this including the 
NMSSM [17, 18], the MNSSM (Minimal Non-minimal Supersymmetric Standard Model) [19] 
and larger models [20] with interesting implications for dark matter. [21] For concreteness and 
the availability of dominant 1-loop and 2-loop corrections to the higgs sector via the code 
NMHDECAY, [22] we choose to study the NMSSM. 

Adding a Higgs singlet is attractive for several reasons. Most interesting, perhaps, is that it 
provides an elegant solution to the /i-problem present in the MSSM [23]. Additionally, the "little 
fine tuning problem" , which results in the MSSM from the lack of a detection of a CP-even Higgs 
at LEP II, is less severe within the NMSSM [7, 24], and is completely absent if the lightest CP- 
odd Higgs is light enough to allow H AiAi decays [7]. Thirdly, baryogenesis considerations 
leave the MSSM in disfavor, requiring the right handed stop squark to be lighter than the top 
quark and the Higgs lighter than about 117 GeV [25]. Recent studies of baryogenesis within 
the NMSSM indicate that parameter points with a light singlet Higgs and a corresponding 
light neutralino are favored [26]. Finally, the domain wall problem [27] in the NMSSM can be 
avoided by the introduction of appropriate non-renormalizable Planck-suppressed operators, 
and imposing a discrete R-symmetry on them. [28] 

In the NMSSM, the physical CP-even and CP-odd Higgs states are mixtures of MSSM-like 
Higgses and singlets. The lightest neutralino therefore has, in addition to the four MSSM com- 
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ponents, a singlino component which is the superpartner of the singlet Higgs. The eigenvector 
of the hghtest neutrahno, Xi? in terms of gauge eigenstates is: 

X? = e„i^° + e,^° + ewW^ + e^^ + e,^, (1) 

where eu, are the up-type and down-type higgsino components, ew, are the wino and bino 
components and is the singlet component of the lightest neutralino. 
Likewise, for the lightest CP-even Higgs state we can define: 



Hi 




(2) 



Here, 3f? denotes the real component of the respective state, and we take vacuum expectation 
values to be those of the complex states (e.g. v = \Jv^l^^^i^ ^ 174 GeV). 
Lastly, we can write the lightest CP-odd Higgs as: 

Ai = cos OaAmssm + sin OaAs, (3) 

where As is the CP-odd piece of the singlet and Amssm is the state that would be the MSSM 
pseudoscalar Higgs if the singlet were not present. 6 a is the mixing angle between these two 
states. There is also a third imaginary linear combination of H^, and S that we have 
removed by a rotation in /5. This field becomes the longitudinal component of the Z after 
electroweak symmetry is broken. 

The NMSSM can contain either an approximate global U{1) R-symmetry in the limit that 
the Higgs-sector trilinear soft SUSY breaking terms are small, or a f/(l) Peccei-Quinn symmetry 
in the limit that the cubic singlet term in the superpotential vanishes [29]. In either case, one 
ends up with the lightest CP-odd Higgs boson, Ai, as the pseudo-goldstone boson of this broken 
symmetry, which can be very light. In some regions of the NMSSM parameter space, one can 
also get the lightest CP-even state. Hi, to be very light as well. This is discussed in more detail 
in Sec. IV. As shown in Sec. HI, it is easy to get a light largely singlino LSP in the U{1)pq 
symmetry limit. 

While we confine our analysis to the NMSSM, it should be noted that such symmetries are 
generically present in other singlet models such as the Minimal Non-minimal Supersymmetric 
Standard Model (MNSSM) [19]. The combination of a light Ai and a light neutralino is 
not uncommon in a wide class of models with extra singlets and/or extra gauge groups [20]. 
Implications of such models for the relic neutralino density have been considered in [30]. 

The NMSSM is defined by the superpotential 

XH^H.S+'^S'' (4) 



4 



^ Singlino LSP 




FIG. 1: On the left, we show regions of X-n parameter space for which the Xi is singhno-hke (defined 
by > 0.5) and bino-hke (defined by < 0.5). On the right, we plot niA vs. m~o for singlino-like 
neutralinos with > 0.9. Each point shown is consistent with all LEP constraints. 



and associated soft-supersymmetry-breaking terms 

XA^H^H.S + ^A^S\ (5) 

where the hatted objects are chiral superfields and unhatted objects are their scalar components. 
An effective /i parameter (as defined by the superpotential form fiHuHd of the MSSM) is 
generated from the first term of Eq. (4) when (S) = x is non-zero: fi = Xx. We follow the sign 
conventions for NMSSM parameters of Refs. [31, 32] in which A and tan/5 = v^/vd are positive 
while K, A\ and can have either sign. 



III. LIGHT NEUTRALINOS IN THE NMSSM 



In the basis x° = (^^-^i? —^-^25 V^^? V^s); the tree-level neutralino mass matrix takes the form 



1 


Ml 










M2 












V2 




gr"d 


a2Vd 



















-/i 



91-" d 

g2Vd 






-Xvu 



(6) 



In the above, the upper 4x4 matrix corresponds to Ai-^f^^. From the lower 3x3 matrix, we 
find that if Xvu^d = {fJ'/x)vu,d are small compared to and/or 2\k,x\ then the singlino decouples 
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from the MSSM and has mass 



"Isinglino — VA^f^ + 4k^x2 = Y li^V^ jx^ + A^K^X^ , (7) 

as found from [M~o]55. Thus, if 2|kx| and \v are both < Mi, M2, then the hghtest neutrahno 
will tend to be singlino-like [33]. Since |x| is typically substantial (given that A < 1 and 
= \\x\ must be substantial to satisfy chargino mass limits), a singlino-like Xi (formally 
defined by < 0.5) emerges mainly for small k. In fact, for very small A, |x| must be quite 
large and thus the singlino will be the LSP only if \k\ is also very small; otherwise 2\kx\ 
would exceed one or more of the typically moderate Mi, M2, values considered here and the 
singlino would not be the LSP. For larger A (e.g. > 0.3), |x| need not be extremely large and 
the singlino LSP condition 2\kx\ < Mi,M2, |/i| can hold for slightly larger These behaviors 
can be seen in Fig. 1 obtained by scanning using NMHDECAY 1.1 [22]. NMHDECAY tests 
for theoretical consistency of the model and for consistency with LEP constraints on the Higgs 
sector, neutralinos and the chargino). It also includes radiative corrections to the tree- level mass 
matrix that are often quite important for small As expected, the neutralino is singlino- 
like (e^ > 0.5) when \k,\ is small. Consistent solutions are found primarily in two regions of 
parameter space — one at small A with very small \k\, and another at large A with slightly 
larger \k\ allowed, see Fig. 1. Since tan/3 also induces singlino mixing, the singlino points at 
large A also have small tan (3 < 4 while the points at small A can have any value of tan (3. 

For |k| not close to zero, bino-like Xi's can easily emerge for small values of Mi. In this case, 
the bino does not have a large degree of mixing with the other neutralinos and the LSP mass 
is nearly fixed to Mi. 

We will find that a light neutralino which is mostly bino or a combination of bino and 
singlino, with a small admixture of higgsino, can generate the observed dark matter density 
and evade all relevant collider constraints. 



IV. LIGHT CP-ODD HIGGS BOSONS IN THE NMSSM 



After removing the CP-odd degree of freedom that is absorbed in giving the Z its mass, the 
remaining CP-odd states have the squared-mass matrix 

^2 ^(^,{A, + Kx) Xv{A,-2kx) 

^ \\v{Ax-2kx) {2\k + ^)v'^ sm2(3 - ?,xkA^ 

where = i;^ + wj. For physically acceptable solutions, the lightest state must have m\_^ > 
0. In addition, the lightest CP-even Higgs boson and the charged Higgs boson must have 
positive mass-squared. To avoid spontaneous CP-violation several other conditions must be 
satisfied [31]. In our conventions these are as follows. 
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• For K > 0, we must have one of three situations: 

1. sign(/i) = sign(AA) = -sign(A«,); 

2. sign(/i) = -sign(AA) = -sign{A^,) with \Af,\ > 3XvuVd\Ax\/ {-\xAx\ + kx^), where 
the denominator has to be positive; 

3. sign(/i) = sign(AA) = sign(yl«;) with \A^\ < 3XvuVd\Ax\/{\xAx\ + kx^). 

• For ft < 0, a CP-conserving minimum requires 

1. sign(;u) = sign(AA) = sign(A^) with \A^\ > 3\vuVd\Ax\/{\xAx\ - kx'^). 

To find a model which has a hght CP-odd Higgs boson, we can require that one of the 
U{l)ii or U{l)pQ symmetries approximately holds. The U{l)ii symmetry appears in the limit 
that the trilinear terms and Ax vanish. This is well motivated from models of gaugino- 
mediated SUSY breaking [34] in which trilinear terms are generated radiatively and, therefore, 
are suppressed relative to the gaugino masses by a loop factor of An. One would expect Ax 
to be smaller than the gaugino masses by a factor of An and to be smaller by a factor 
of IGtt^ because S is not charged under gauge symmetries, and only receives a trilinear term 
at two-loops. The small trilinear terms are also radiatively protected and remain small when 
evolved via RGEs from the SUSY breaking scale to the weak scale [29]. In this limit, the 
lightest CP-odd Higgs is a pseudo-goldstone boson of the broken f/(l)_R symmetry and has a 
mass of m\_^ ~ —SkA^x in the large tan (3 or large |x| limits. Alternatively, we can make the 
substitution, x = fi/X, and write this as ~ —SjA^fi. 

More generally, in the limit of small Ax and A^ one finds 

_ X '^^ sin^ 2/5 

tan^A ~ . . cos Oa ^ ■ 2^/? , — 2 ' 

wsm2p v'sm 2p + x^ 



and 

~ ^ 

^2 _|_ y2 gjj^ 2/3 



, ^XAxv'^xsm23 — SkAi^x^ , , 



Since |x| > wis preferred, | cos Qa\ is typically small at small to moderate tan /?, with cos^ 6*^ ^ 
at large tan/5. If we only take — ^ 0, one finds the results 

" 1 +AA/(ra) " 



X 

tan^A 

V sm 2/5 

and 



1 - Ax/{2kx) 



'11 



2 fAAAi'^xsin 2/5 

""^1 ~ + sin2 2/5 + Aax/^ ' ^ ' 

valid whenever the numerator of the preceding equation is much smaller than the square of 
the denominator, as for example if tan/5 — > 00 or x is large. Again, cos6'a will be quite small 
typically and the Ai relatively singlet like. In practice, this limit is very frequently applicable. 
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The U{l)pQ symmetry appears in the hmit that k, vanishes (and therefore the soft SUSY 
breaking term also vanishes) and also results in a light Ai. To leading order, one finds 

, n 2x 2 sin^ 2/3 

tan^A , cos Oa ~ 2 ■ , a 2 (^^j 



and 



2 6KX^{3\v'^sm2p -2A^x) 



"^^^ Ax^ + v^sin^2(3 " ^'^^ 

For |x| > V, I cos^^l is small for moderate tan/3 and approaches at large tan/3. 

It is useful to note that if A is small (implying large \x\), then singlet mixing in both 
and Ai\ is small. If 2\kx\ < Mi,M2, \fi\ and 2fi{Ax + kx) / sin 2/3 > — SkA^x then the Xi ^i-nd 
Ai will both be singlet in nature. In particular, for small \k,\ both the Ai and the Xi can easily 
be singlet-like. At large A, the Ai can have a more mixed nature (cos6'yi tends to be larger) but 
as we have seen a moderately-singlino Xi is still allowed despite the somewhat larger mixing in 
the neutralino mass matrix. 

In either of the Ax, A,^ ^ oi k ^ cases, a light Ai is technically natural since it is 
protected by an approximate symmetry. From an effective field theory perspective, the small 
terms that break the symmetry will not receive large radiative corrections. It is technically 
natural for the Ai, Hi or Xi to be very light as a result of U{1)r and/or U{1)pq symmetries. 

The fermion Yukawas break the U{l)fi symmetry, leading to contributions arising at one- 
loop for the Hu and components of the Higgs sector, and at two-loops for the S component. 
However the radiative corrections to the singlet component are proportional to either A or k, 
and thus are suppressed for small values of A, k. These symmetries therefore result in the 
hierarchy rrih ^ Ax ^ A,^. 

It will be helpful in understanding dark matter relic density issues to examine whether or 
not a light (singlet) Ai can decay to a pair of nearly pure Xi singlinos. From the right hand 
plot of Fig. 1, we observe that it is impossible to obtain m^^ > 2m~o when the LSP is nearly 
purely singlino (e^ > 0.9). Using the fact that a highly singlino Xi is achieved by taking A to 
be very small (so as to remove mixing in M~^), implying large we are able to analytically 
understand this in two cases: (i) \kx\ moderate in size and (ii) \k,x\ small. For moderate \kx\, 
the inability to satisfy the mass requirements for the decay Ai XiXi stems from the inability 
to simultaneously satisfy m|g > and ■m\_^/4m'~o > 1. Here, S3 is the third CP-even (largely 
singlet) Higgs mass eigenstate as defined in [17] and is the lightest CP-even Higgs state in the 
limit of interest. For small A and finite \nx\, we have 

m~o ~ 2kx , m\_^ ~ — SkA^x , (15) 
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the latter requiring kA,^ < 0. Further, if one expands the CP-even mass matrix in the large |x| 
limit, holding fi and kx fixed, one finds [17] (after correcting for differences in sign conventions) 



ml^ = iK^x^ + kA^x + ^ ' ' 



2.,2 r,, 1 / /I. \ n2 ^2^,2 



, 2fi; + — )sin2/3 
X 2 V X 



(-2k+ cos"^ 2[3 
V -1 ^ (16) 



4^2x2 - 2me^ 

sm 2/5 



where As = A\ — nx. For fixed with |x| very large, the last two terms approach zero and we 
have m|g ~ An'^x'^ + nA^^x which is positive only if —kA,^ < Ati^x. However, in the same limit, 
the mass condition for the Ai XiXi? written as Tn\^ > 4m~o, becomes —kA,^ > (16/3)fi:2x. 
These two conditions cannot be simultaneously satisfied, and thus the decay Ai XiXi is 
not allowed for a pure singlino in the large fixed kx limit. Some admixture of bino, and 
therefore moderate A and \x\ must be required for this decay to be open. 

The |k| — > case (the Peccei-Quinn symmetry limit) at small A is defined by 

I I ^ A I^aI \A^\ \fi\ V \ 

|k| < C A, , ,1—^,1—^], v<t:\x\. (17) 

\ V V \x\ \x\ J 

In this limit, Eq. (7) implies m~o ^ A^f^ = fi^v'^/x'^. Meanwhile, for \k,x\ — > and |x| large 
(A small), it is easily seen that m\_^ oc l/|xp. Thus, once again, the Ai —>■ XiXi decay is 
disallowed. 

The fact that ttlai < 2m~o in these singlino limits implies that a singlino Xi is disfavored 
cosmologically. This is because m^^ ~ ^mjo is required to enhance the annihilation cross 
section to the level needed to obtain the correct relic density. Therefore some bino mixing is 
required to get an appropriate relic density. This, in turn, requires Mi to be small as well. 

Finally, we should note that the constraints on a light Ai are rather weak. This is because 
direct searches at LEP [35] require a light Ai to be radiated off of a quark or a tau lepton 
and, due to the small fermion Yukawa couplings, bounds are only obtainable when the Ai 
coupling to fermions is enhanced by tan/5. However, in the U{1)pq and U{1)r symmetry limits 
discussed earlier, cos 9a (the non-singlet part of Ai) is proportional to sin 2/3 so that the product 
tan j3 cos 6a remains modest in size. If a light Ai exists and is near in mass to the t] (547 MeV), 
it may be discovered via invisible decays of the t] at low energy lepton colliders. [36]. This mass 
range is extremely interesting if a light Ai exchange is the explanation for the recent galactic 
511 keV line from the INTEGRAL/SPI experiment [10]. 

We now give some additional remarks concerning the singlet and non-singlet Ai possibilities. 

A. Models with a Singlet-Like Ai 

As we have seen, a singlet-like Ai [cosOa ^ 1) is extremely easy to obtain by making some 
combination of IkI, \AJ and \Ax\ small. Indeed, the mass of the Ai can be driven to zero at 
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tree level. Radiative corrections increase this mass, however. The dominant source of these 
radiative corrections to the singlet mass is from Standard Model couplings since the number 
of degrees of freedom is much larger in the MSSM than in the singlet supermultiplet. These 
radiative corrections are therefore proportional to A, since the A superpotential term is the only 
coupling connecting the singlet with the rest of the MSSM. Therefore, if the light singlet mass 
is to be radiatively stable, A must be small. A being small also has the effect of reducing the 
mixing with the singlet component in both the CP-even mass matrix and the neutralino mass 
matrix. All terms which mix the singlino with the MSSM neutralinos and the singlet S with Hu 
and Hd are proportional to A. We find A < 0.1 to be natural, with larger values of A requiring 
an increasing amount of cancellation between the various radiative contributions to its mass. A 
being this small necessarily implies that the singlet vacuum expectation value, |a;|, is large since 
H = Ax. Chargino searches generally imply > 100 GeV, leading to |x| > 1 TeV for A < 0.1. 
Furthermore, with all four of A, k, A,^ and Ax small in magnitude, the entire supermultiplet is 
light and Ai, Hi (largely the singlet-like S3) and Xi tend to be nearly degenerate. 

B. Models with an MSSM-Like Ai 

An MSSM-like (non-singlet) Ai (cos 6'^ ~ 1) can also be obtained, but is subject to more 
stringent constraints. If cos 6'^ ^ 1, couplings of the Ai to down-type fermions go like 
cos 0^ tan/?, therefore phenomenological constraints become significant at large tan/5. If such 
an Ai is very light, it will be further constrained by rare decays such a.s K ^ iruu and T 7X 
as discussed in Sec. VE. 

V. CONSTRAINTS 

In this section, we will consider a series of constraints which may be relevant to light neu- 
tralinos and/or a light CP-odd Higgs bosons in the NMSSM. Except for the LEP and T decay 
limits, most of the constraints discussed below are easily avoided by appropriate choices of 
SUSY parameters to which our dark matter calculations are not sensitive. 

A. LEP Limits 

If the lightest neutralino is lighter than Z decays to neutralino pairs may violate the 

bounds obtained at LEP for the Z's invisible decay width. In particular, we require P^^-o-o < 
4.2 MeV, which corresponds to one standard deviation from the measured neutrino contribution. 
Since binos, winos and singlinos do not couple to the Z, this constraint can only limit the 
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higgsino components of the lightest neutrahno. In the mass range we are most interested in 
here {rn~o < 20 GeV), this constraint is satisfied for all models with — < 6%. 

Direct chargino searches also limit the wino component of the lightest neutralino. This is 
due to the fact that if the lightest neutralino has a significant wino component, then it will 
have a mass that is a significant fraction of the chargino mass (with near degeneracy if the Xi 
is mainly wino). 

In combination, these constraints imply that a very light Xi must be dominantly a linear 
combination of bino and singlino. 



B. The Magnetic Moment of the Muon 



The one-loop contribution to the magnetic moment of the muon from a light neutralino 
comes from a triangle diagram with a smuon along two sides and the neutralino around the 
third. This contribution is given by [37]: 

:ro m,, 



^<^ = i^ E Y23^(i-^r+i-^r)i^f(-j+^Re[«]i^2"(-. 



-L"" m=l,2 



where = V2gieBXm2 + y^e^X^i, = {g2ewX^2 + gi€B)X*^^/ y/2 - y^euX*^^^ = 
g2m^/ {\/2'm\Y cos (3) and X^,™ are elements of the unitary matrix which diagonalizes the smuon 
mass matrix. The functions, F^, are defined by: 

Fi{xm) = j—^ — —(l-6xm + 5xl^ + 2xf^~6xl^\nx^\ (19) 

\ J- Xjji ) \ / 

F2{^m) = —^^(l-xl + 2x^\nxm), (20) 
where Xm = m~o/mf, . For a bino-like neutralino, this reduces to: 



,2 



(21) 



where is the real part of the off diagonal elements of the smuon mass matrix, 
3?[m^(A^ — jj* tan/3)]. This further simplifies to: 



2 3 X / 200GeV \Y ^tan/3-A^ 
,10 GeV A mr, ) \ 1000 GeV 



(5ay ^2.3 X 10-^M^^1( 1 ( ' n . (221 

In addition to this contribution from a light neutralino, a light CP-odd Higgs can contribute 
non- negligibly to 5a^ through both one-loop and two-loop processes [38]. The one-loop con- 
tribution (corresponding to a triangle diagram with a muon along two sides and the CP-odd 
Higgs along the third side) is given by: 



gl ml cos^ 6a tan^ P La 
32 TT^ 



^ A Hoop _ ^2 "> ^""^ — A" , . 
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where the function, La, is given by: 

= /' -^'^^ .24) 

■m\ Jo ("^^/''^a) + (1 — a^) 

Numerically, this function yields La = -0.032, -0.00082 and -0.000014 for 1, 10 and 100 GeV 
CP-odd Higgs bosons, respectively. Thus we arrive at: 

^^AUoop ^ LACOs^eAtan^ P x 2.7 x 10-^ (25) 

which is a negative contribution due to the sign of La- 

The contribution from two-loop diagrams involving a heavy fermion loop is given by: 

.,A2ioop _ 9lmlacfQ)x)Lf 

^^^^ ~ 32m2^7r2 ' ^^""^ 

where C/ is the color factor of the fermion in the loop (3 for quarks, 1 for leptons), Qf is 
the electric charge of the fermion, a ~ 1/137, Xf = cos 6*^ tan/? for up-type fermions and 
cos cot /3 for down- type fermions and Lf is given by: 

^ 2m\ lo X (1 — x) — (m^/m^) ( rnj/m\ ) 

For a top quark loop, this function yields the values Lt = 6.2, 3.9 and 1.7 for 1, 10 and 100 
GeV CP-odd Higgs bosons, respectively. For a bottom quark loop, these values are L^ = 2.5, 
0.59 and 0.038. Numerically, these contributions are: 

2 loop _ ^Qg2 ^^^2 ^ X 2.6 X 10"^^ + Lb cos^ 9a tan^ x 6.6 x 10~^l (28) 

Combining the results of the one and two-loop contributions from the light Ai, we arrive at the 
following: 

1+2 loop ~ _7 X 10^" X cos^^Atan^/? for rriA, = 1 GeV, 
^^Ai+2ioop ^ ^ 7 ^ ^Q-12 ^ cos^^Atan^/? for rriA, = 10 GeV. (29) 

It is somewhat difficult to know how best to interpret the current status of the measurement 
of the muon's magnetic moment. Using e~^e~ data, the measured value exceeds the theoret- 
ical prediction by Sa^{e~^e~) = [23.9 ± 7.2had-io ± 3.5ibi ± 6cxp] x 10~^°, where the error bars 
correspond to theoretical uncertainties in the leading order hadronic and the hadronic light-by- 
light contributions as well as from experimental contributions. Combined, this result is 2.4cr 
above the Standard Model prediction. Experiments using t~^t~ data, on the other hand, find 
5a^(r"'"r~) = [7.6 ± 5.8had-io ± 3.5ibi ± 6cxp] x 10""^°, which is only 0.9a above the Standard 
Model prediction [39]. 
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Comparing the expression shown in Eq. (22) for a hght neutrahno to these experimental 
results illustrates that only in extreme models, with a combination of small m^, large tan/3 
and large /i is there any danger of exceeding these bounds with a light neutralino. For fairly 
moderate choices of parameters, i.e. ~ 200 GeV, tan/3 ~ 20 and /i ~ 500 GeV, the 
experimental values can be matched. Since is not currently known and the dark matter 
scenarios we consider are not sensitive to m^, we will not consider this constraint in our dark 
matter calculations. 

The contribution from a light CP-odd Higgs also should not violate the 5a ^ constraint. 
If one considered the Ai contribution alone, one might conclude that cos tan /3 is strongly 
limited in the case of small m^^ . However, contributions to 5a^ from other sources such as the 
charged Higgs, charginos, and sfermions can easily overwhelm or cancel any contribution from 
a light Ai. Furthermore, LEP and other indirect limits such as T decays (discussed in Sec.VE) 
constrain cos tan /3 to be small, so it is generally not possible to see a large enhancement in 
5a^. Finally we note that if the Ai and Xi ai'e both light, as considered here, their contributions 
to 5a^ are of opposite sign, and can cancel. 

Thus, we do not explicitly include the constraints in our computations. Their inclusion 
would only become appropriate if a specific model for soft-SUSY-breaking is being considered. 

C. Rare Kaon Decays 

The Tx^vv branching ratio was recently measured by the E787 and E949 experiments 

to be BR(K^ -n^uu) = (1.47tJ;|g) x 10~^°, which is nearly twice the value predicted in the 
Standard Model, (0.67t[J:^?) x 10"^° [40]. A CP -even Higgs boson lighter than a few hundred 
MeV can contribute to this branching ratio via a triangle diagram involving bosons on 
two sides, and an up or charm quark on the third. This contribution is suppressed by C,u and 
^d- CP-odd Higgs bosons, on the other hand, cannot contribute to this process at the one 
loop level since the vertex involving W^s and the Ai is W^W^AiAi and, therefore, the leading 
contribution to 7r++ invisible has four x^'s in the final state. This requires a x? lighter 

than 88.5 MeV, which is lighter than the range we consider in this study. 

Other rare kaon decays such as —>■ e^e~ and 7r+e+e~ are similarly unconstraining 

for a light Ai, but potentially important for a light Hi for the same reasons. 

A recent study [41] analyzed this in detail and concluded that extremely light m^^ < 2m^ 
can be ruled out. However, this can be evaded if |k| is small enough. 
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D. Rare i?-Meson Decays 



The transitions h ^ and h are usually considered sensitive probes of supersym- 

metry, however both are flavor changing, while a light Ai and Xi are not flavor changing by 
themselves. These and other flavor changing processes involving a light Xi propagator can al- 
ways be suppressed by making the appropriate squark or slepton mass heavy since the relevant 
diagrams must involve a ffxi vertex. Processes involving a light Ai may be suppressed at 
one- loop by assuming the Minimal Flavor Violation mechanism [41, 42]. A recent study of the 
B meson decays b — * 57, b sAi, and b sl^l~ in the NMSSM concluded that Ai masses 
down to 2me cannot be excluded from these constraints [41]. 

Another rare 5-decay is K^vv. This process also necessarily involves a quark flavor- 

changing vertex. A diagram in which the light Ai couples to the must involve two 
Ais and two W^^s unless CP is violated, severely limiting the set of processes to which it can 
contribute. Diagrams where the light Ai couples to the fermion also must have a to change 
the quark flavor and also receive a factor of cos^a at each fjAi vertex, strongly suppressing 
the Ai contribution for the scenarios we focus on, all of which have small cosOa- 

E. Upsilon and J/^ Decays 

The vector resonances J /'^ and T may decay radiatively into an Ai and a photon if Ai is 
sufficiently light. There are two experimental limits on this process: firstly when the Ai decays 
invisibly or is long-lived enough to leave the detector volume [43], and secondly when the Ai 
decays to Standard Model particles [44]. This width, relative to the width to muons at leading 
order is given by [45]: 



where V is either J/\I' or T and X = cos6'Atan/3 for T and X = cos 6 a cot (3 for J/\E'. The Ai 
is often referred to as the axion in this literature. Eq. (30) is also applicable for a light CP-even 
Hi, with X = id/ cos (3 for y = T and X = ij sin/3 for V = J/^. 

It is usually claimed that a light MSSM A is ruled out if it is light enough so that both the T 
and J/"^ can decay to it {niA ^ 3.1 GeV), due to the observation that r(J/\E' — > 7^4) x r(T 
7^4) is independent of tan (3. However, within the NMSSM, this product is proportional to 
cos'^ 6a, which may be small. 

The best existing measurement of T ^ invisible+'-f is from CLEO [43] in 1995. Significantly 
more data has been collected on the T(IS') resonance that could be used to improve this 
measurement, however. Modern i?-factories such as BaBar and Belle can also produce the 
T(IS') in Initial State Radiation to improve this measurement. 




(30) 
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FIG. 2: The branching ratio for TflS") — > 7X1X1 via 3-body decay (i.e. either rriA < 2m~o or 
ruAi > nir) is plotted vs. the LSP mass (left) and relic density flh'^ (right). All points shown are 
consistent with all LEP constraints. Points marked by an x are excluded by one of: T — > 7X1X1 
(3-body decay) (that which is plotted); T — > 7^1 (2-body decay) with Ai XiXi (2-body decay); or 
T — > 7A1 (2-body decay) where the Ai decays visibly. 

In Fig. 2, we illustrate the correlations between BR{T 7XiXi) via 3-body decay (i.e. 
not T — > 7A1 with Ai invisible or visible), m~o and the relic density Qh"^ (the calculation 
of which is discussed in the following section). Only higgs exchange is included in these relic 
density values. Sub-leading Z and sfermion exchanges would further decrease the relic density of 
points with very large Qh'^. The left-hand plot shows that a significant fraction of the parameter 
choices such that T —* 7X1X1 is allowed are eliminated by the experimental constraint on this 
mode, with additional ones being eliminated by the constraints on the 2-body T '-/Ai decay 
mode. But, many are not excluded, especially those with a bino-like Xi- Improvement in the 
experimental sensitivities to BR{T —>■ 7X1X1) ^"^^ BR{T 7^1) will further constrain the 
light Xi scenarios considered here, or could yield a signal. The right-hand plot of Fig. 2 shows 
that there are many parameter choices that yield BR{T 7XiXi) ^"^^ BR{T 7^1) below 
the experimental limits while simultaneously predicting a relic density roughly consistent with 
observation. We observe that this dual consistency can achieved for either a bino-like or a 
singlino-like lightest neutralino. 

The points on the left side of the right frame of Fig. 2 undergo T 7XiXi dominantly via 
a CP-even, mostly-singlet scalar. Hi, which mediates this interaction. When the Ai becomes 
light and mostly singlet, it often brings the CP-even scalar and singlino down in mass as well. 
For these points, the two body decay, T 7^1, followed by the decay, Ai XiXi? is also just 
below the experimental limit. 
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VI. ANNIHILATION CROSS SECTION AND RELIC ABUNDANCE 



The calculation of the neutralino annihilation cross section and relic abundance in the 
NMSSM is only slightly modified from the case of the MSSM. First, the diagonalization of 
the 5x5 neutralino mass matrix of the NMSSM yields different LSP compositions for given 
choices of input parameters (Mi, /i, etc.). Secondly, annihilations can occur through the ex- 
change of a Higgs boson with a significant singlet component. On one hand, this weakens the 
respective couplings. On the other hand, much lighter Higgses can be considered, as collider 
constraints are weakened. 

For the range of masses we are considering, the only final states available for the annihilations 
of light neutralinos are fermion pairs. This process can occur through s-channel Higgs exchange 
(both CP-even and CP-odd), s-channel Z exchange or t and w-channel sfermion exchange. 
With LEP constraints limiting sfermion masses to vtij > 100 GeV, neutralinos lighter than 
approximately 25 GeV cannot annihilate efficiently enough through sfermions to yield the 
measured relic density. Similarly, Z exchange cannot dominate the annihilation cross section 
for light neutralinos. Therefore, we focus on the process of Higgs exchange. 

The squared amplitudes for the processes, XiXi ^ ^ ^ f f and XiXi H ^ f f, averaged 
over the final state angle are given by [46]: 




~ {s-mlY + mlTl 167r V s' ^ ' 

where the labels A and H denote a CP-odd and CP-even Higgs, respectively. Here, C^j^, ^//h' 
C~o~o^ and ^^0^0^ 

are the fermion-fermion- Higgs couplings and the neutralino-neutralino-Higgs 
couplings, and rriA H and Ta,h are the Higgs masses and widths. In the NMSSM case, we will 
be considering only A = Ai and H = Hi, the lightest of the CP-odd and CP-even states, 
respectively. The relevant couplings are then given by: 

C~o~o^ = cos9a {g2ew - gi(^B){(^dCosP - euSinP) + VlXesieuSinP + eacosP) (33) 
+ sin 6^^1/2 Ae„erf - Ke^ 



C 



ffA 



cos6'^ tan/3, f = d,s,b,l 
^cos6'ACot/3, f = u,c 



(35) 
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c 



ffH 



y/2v COS P 



, / 
/ 



(i, s, &, / 



(36) 



V2d sin fi ' 

We expect eV-MeV \i A = is mostly singlet and 1-10 MeV otherwise. 

Similarly, we expect Yu ^10 eV-100 keV if H\ is mostly singlet and Yu ^ keV-MeV ii H = Hi 
is mostly non-singlet. These widths are strongly affected by the many kinematic thresholds 
due to hadronic resonances with masses less than 10 GeV. Therefore, any computation of the 
relic density is inherently limited by our ability to compute hadronic form factors and sum over 
hadronic decays which may be on-shell and may enhance the annihilation. We require only 
that the relic density is (9(0.1). There is sufficient parameter space to make the relic density 
precisely the value measured by WMAP when all hadronic corrections are taken into account. 
In our computations, we neglect the widths since they are very small compared to the masses 
considered. Of course, one could always tune 2m~o to some hadronic resonance or threshold in 
order to drastically increase the cross section and thus reduce the thermal relic density, but we 
do not employ such precision tuning. 

The squared amplitudes of Eqs. (31) and (32) can be used to obtain the thermally averaged 
annihilation cross section [47]. Using the notation sq 



4:771^0, we have 

Xi 



{av) 



mi 



Xi 



3T 



Uj{So) 
Xi 



2uj'{so) 



T + 0{T^ 



(37) 



s— >4mi +6m~nj 
^0 y. 



where T is the temperature. Keeping terms to zeroth and first order in T should be sufficient 
for the relic abundance calculation. Writing this as an expansion in x = T/m~o, {av) = 
a + bx + (9(a:^), we arrive at: 



a 



glcfrnj cos^ Oa tan2 f] "^|o ^l-mj/m^o 



(4m2 



m 



2 \2 



+ m\Y\ 



X 



eu(eiy — efitan^iy) sin/3 + ed(evi/ — e_Btan6'iy) cos/3 



-|- — es(e^t sin/3 + cos/3) 

92 



tan 9 a 

92 



V2(Ae,. 



(38) 



^XX^A^ff 



0, 



(39) 



where Cf is a color factor, equal to 3 for quarks and 1 otherwise. For this result, we have 
assumed that the final state fermions are down-type. If they are instead up-type fermions, the 
couplings used must be modified as described above. 
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FIG. 3: The CP-odd Higgs mass required to obtain the measured reUc density for a light neutrahno 
in the MSSM. Models above the curves produce more dark matter than in observed. These results are 
for the case of a bino-like neutralino with a small higgsino admixture (e^ = 0.94, = 0.06). Results 
for two values of tan/3 (10 and 50) are shown. The horizontal dashed line represents the lower limit on 
the CP-odd Higgs mass in the MSSM from collider constraints. To avoid overproducing dark matter, 
the neutralino must be heavier than about 8 (22) GeV for tan/3 = 50 (10). 

We have not written the result for CP-even Higgs exchange because the low velocity term in 
the expansion is zero: a^^_,ff_,j-f — 0. Although the 6-terms can, in principle, contribute to the 
freeze-out calculation, in the computations here such contributions do not have a significant 
impact. 

The annihilation cross section can now be used to calculate the thermal relic abundance 
present today. 



where is the number of relativistic degrees of freedom available at freeze-out and xpo is given 



For the range of cross sections and masses we are interested in, xpo ~ 20. 

As a benchmark for comparison, we consider a light bino which annihilates through the 
exchange of an MSSM-like CP-odd Higgs (cos^a = 1)- The results for this case are shown 
in Fig. 3. In this figure, the thermal relic density of LSP neutralinos exceeds the measured 
value for CP-odd Higgses above the solid and dashed curves, for values of tan /5 of 50 and 10, 




(40) 



by: 




(41) 
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FIG. 4: We display contours in niAi ~ rn~o parameter space for which Eq. (40) yields O/i^ = 0.1. 
Points above or below each pair of curves produce more dark matter than is observed; inside each set 
of curves less dark matter is produced than is observed. These results are for a bino-like neutralino 
with a small higgsino admixture (e^ = 0.94, = 0.06). Three values of tan/3 (50, 15 and 3) have 
been used, shown as solid black, dashed red, and dot-dashed blue lines, respectively. The dotted line 
is the contour corresponding to 2m~o = tjia- For each set of lines, we have set cos^ 6a = 0.6. The 
tan /? = 50 case is highly constrained for very light neutralinos, and is primarily shown for comparison 
with the MSSM case. 

respectively. Shown as a horizontal dashed line is the lower limit on the the MSSM CP-odd 
Higgs mass from collider constraints. This figure demonstrates that even in the case of very 
large tan (3, the lightest neutralino must be heavier than about 7 GeV. For moderate values of 
tan (3, the neutralino must be heavier than about 20 GeV. 

In Fig. 4, we show how this conclusion is modified within the framework of the NMSSM. Here, 
we have considered a CP-odd Higgs which is a mixture of MSSM-like and singlet components 
specified by cos^ 9a = 0.6 and a neutralino with composition specified by = 0.94 and 

= 0.06. These specific values are representative of those that can be achieved for various 

NMSSM parameter choices satisfying all constraints. For each pair of contours (solid black, 

dashed red, and dot-dashed blue), the region between the lines is the space in which the 

neutralino's relic density does not exceed the measured density. The solid black, dashed red, 

and dot-dashed blue lines correspond to tan/9=50, 15 and 3, respectively. Also shown as a 

dotted line is the contour corresponding to the resonance condition, 2m~o = rriA- 

xi 

For the tan/5=50 or 15 cases, neutralino dark matter can avoid being overproduced for any 
Ai mass below ~ 20 — 60 GeV, as long as m~o > m?,. For smaller values of tan 3, a lower limit 
on niAi can apply as well. 

For neutralinos lighter than the mass of the 6-quark, annihilation is generally less efficient. 
This region is shown in detail in the right frame of Fig. 4. In this funnel region, annihilations 
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to cc, T~^T~ and ss all contribute significantly. Despite the much smaller mass of the strange 
quark, its couplings are enhanced by a factor proportional to ta.n/3 (as with bottom quarks) 
and thus can play an important role in this mass range. In this mass range, constraints from 
Upsilon and J/ip decays can be very important, often requiring fairly small values of cos^a- 

For annihilations to light quarks, cc, ss, etc., the Higgs couplings to various meson final 
states should be considered, which include effective Higgs-gluon couplings induced through 
quark loops. In our calculations here, we have used the conservative approximation of the 
Higgs-quark-quark couplings alone, even for these light quarks, but with kinematic thresholds 
set by the mass of the lightest meson containing a given type of quark, rather than the quark 
mass itself. This corresponds to thresholds of 9.4 GeV, 1.87 GeV, 498 MeV and 135 MeV for 
bottom, charm, strange and down quarks, respectively. A more detailed treatment, which we 
will not undertake here, would include the proper meson form factors as well as allowing for 
the possibility of virtual meson states. 

Thus far, we have focused on the case of a bino-like LSP. If the LSP is mostly singlino, it 
is also possible to generate the observed relic abundance in the NMSSM. A number of features 
differ for the singlino-like case in contrast to a bino-like LSP, however. First, the ratio m-o/rriA^ 
cannot be arbitrarily small. The relationship between these two masses was shown for singlino- 
like LSPs in Fig. 1. As discussed earlier, and shown in this figure, an LSP mass that is chosen 
to be precisely at the Higgs resonance, m^^ ~ 2m~(), is not possible for this case: rriAi is always 
less than 2m~o by a significant amount. 

Second, in models with a singlino-like LSP, the Ai is generally also singlet-like and the 
product of tan^ j3 and cos^ 6a is typically very small. This limits the ability of a singlino- 
like LSP to generate the observed relic abundance. The last term in Eq. (38) introduces an 
additional tan^ 6a dependence, however, which effectively reduces the impact of cos 6a on the 
annihilation cross section from four powers to two. But, this last term is suppressed when the 
singlet fraction es is large and eu, are small by the factor of k (which is small for a singlino) 
that multiplies e^. Alternatively, the second to last term in Eq. (38) can also be of importance. 
Overall, the inability to compensate the smallness of the coefficients in Eq. (38) by being nearly 
on-pole implies that annihilation is too inefficient for an LSP that is more than 80% singlino. 

In the following section, we give sample cases for which and m~o are light and Qh"^ ~ 
0.11. These are representative of the many different types of scenarios that are possible and 
include a case in which the Xi is largely singlino. 



20 



VII. SAMPLE MODEL POINTS 



In this section we present specific sample model points of the type we propose. These points 
are obtained using NMHDECAY 1.1. [22] 

The first has a singlet-like Hi, which would have escaped detection at LEP due to this singlet 
nature. In addition, the mass of the more SM-like H2 is beyond the LEP reach. It also has a 
sizable BR{T ^ 'j + Ai) which could be discovered by a re-analysis of existing CLEO data. 
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K 
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tan P 
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-0.00594669 


0.11476 


0.26493 


es 
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4.55841e-26 cm^/s 
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TABLE I: Sample model point #1. 



The second point has an MSSM-like Hi, but due to the presence of the light Ai and the large 
A coupling, this MSSM-like Hi decays dominantly to a pair of Ais [BR{Hi AiAi) = 99.6% 
for this point]. Such an Hi would not be easily detected at the LHC. 
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TABLE II: Sample model point #2. 



The third point has a singlino-like Xi ^^11 as a singlet-like Hi. As for point #1, this point 
has a BR{T ^ 7 -|- Ai) that might be excluded by an appropriate re-analysis of existing data. 
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TABLE III: Sample model point #3. 



VIII. ELASTIC SCATTERING OF LIGHT NEUTRALINOS 



The spin-independent elastic scattering cross section of a light neutralino with nuclei is 
generally dominated by the t-channel exchange of a CP-even Higgs boson. The cross section 
for this process is approximately given by: 

« [^^] {^C,-,MN)) (42) 

where the first sum is over the CP-even Higgs states of the NMSSM and mu are their masses. 
The second sum is over the quark types and < N\qq\N > are the matrix elements over the 
atomic nuclear state. Of course, one must be careful to use the correct form of CqgH which 
differs for up-type quarks versus down-type quarks. In the sum over quark species, the strange 
quark contribution dominates with rris < A^|ss|A^ >~ 0.2 GeV. For a bino-like LSP and any 
one H, Eq. (42) reduces to 

^clastic ~ — : sm &w (edU - e«4dj x 

TTmjj \mp + m~o J 

E ^ < iV|gg|iV > + 5: ^ < N\qq\N > Y. 
If the LSP is singlino-like, on the other hand, the appropriate approximation is 

singlino ^GWz ( ^P^X? ^ ^ 2A^ 6^ COS^ g^4^ ^ 2^ 

E ^ < iV|gg1iV > + E ^ < N\q-q\N > 

A f. COS p J—',^ sm ti 

q=d,s,b " q—u,c 

where, in C~o~Ofj, we have dropped the term containing k since it is expected to be very small 
and the term proportional to e^j which is also likely to be very small. 



22 



In assessing the implications of the above, it is useful to note that LEP limits on a Higgs 
boson with rriH < 120 GeV generally imply 



and for a light Xi LEP limits on invisible Z decays roughly imply e^ rf < 0.06. 

The claim of a positive WIMP detection made by the DAMA collaboration is not consistent 
with the limits placed by CDMS and others for a WIMP in the mass range normally considered 
(above a few tens of GeV). Very light WIMPs, however, scatter more efficiently with light target 
nuclei than with heavier nuclei, which can complicate this picture. For a WIMP with a mass 
between about 6 and 9 GeV, it has been shown that the DAMA results can be reconciled with 
the limits of CDMS and other experiments [16].^ This is made possible by the relatively light 
sodium (A=23.0) component of the DAMA experiment compared to germanium (A=72.6) and 
silicon (A=28.1) of CDMS. 

To produce the rate observed by DAMA, a light WIMP would need an elastic scattering cross 
section of 7xl0~^° cm^ to 2x10"'^^ cm^ (0.7— 2 fb). For the case of a bino-like or singlino-like 
neutralino capable of resolving the DAMA discrepancy, the scale of this cross section is: 



One has Fa = 1 for the bino-like case and Fa = 2A^/(5'| tan^ 6'iy) ^ 0.67 x (A/0. 2)^ for the 
singlino-like case. For tan/3 = 50, A = 0.2 and a Higgs mass of 120 GeV, we estimate a 
neutralino-proton elastic scattering cross section on the order of 4 x 10~^^ cm^ (4 x 10^^ fb) 
for either a bino-like or a singlino-like LSP. This value may be of interest to direct detection 
searches such as CDMS, DAMA, Edelweiss, ZEPLIN and CRESST. To account for the DAMA 
data, the cross section would have to be enhanced by a local over-density of dark matter [16]. 

The cross section in Eq. (44) is small unless tan f] is quite large, in which case the scenario 
will run into difficulty with LEP limits unless cos6'yi is quite small. To explain the DAMA 
result, we can instead require tuh to be small. For instance, with m~o = 6 GeV, rriH = 3 GeV, 
and tan/3 = 10, the DAMA result can be reproduced with (Xciastic ~ 4 x 10"^^ cm^ (~ 4 fb), 
without requiring a dark matter wind through our solar system. It would not be unusual for 
a mostly-singlet Hi to be this light if A is small. In this case the singlet decouples from the 
MSSM and the whole singlet supermultiplet is light. 

^ If a tidal stream of dark matter is present in the local halo, WIMP masses over a somewhat wider range can 
reconcile DAMA with CDMS as well. 




(43) 




0.06. 



(44) 
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For a detailed study of direct detection prospects for heavier neutralinos in the NMSSM, see 
Refs. [31, 32]. We find consistency with their resuhs concerning annihilation through H and A 
resonances. 

IX. EXTREMELY LIGHT NEUTRALINOS AND THE OBSERVATION OF 511 KEV 
EMISSION FROM THE GALACTIC BULGE 

If the LSP's mass is even smaller, below ~ 1 GeV, it may still be possible to generate the 
observed relic density. In this mass range, in addition to annihilations to strange quarks {K^ , 
K'^), final state fermions can include muons and even lighter quarks (tt^, vr"). 

There is a Xi mass range in which neutralinos will annihilate mostly to muon pairs. This 
range is < m~o < + m^o/2, or 106 MeV < m~o < 207 MeV. (The upper limit will be 
explained shortly.) This range of parameter space is of special interest within the context of the 
511 keV emission observed from the galactic bulge by the INTEGRAL/SPI experiment. Muons 
produced in neutralino annihilations will quickly decay, generating electrons with energies of 
~ m~o/3, which may be sufficiently small for them to come to rest in the galactic bulge before 
annihilating. 

The upper limit above derives from the fact that the XiXi annihilations should not create 
many vr^'s. In this way, we avoid gamma ray constraints from EGRET. If we assume that 
the annihilation mediator is the CP-odd Ai, XiXi ~^ Ai ^ pions is only possible if 2m~o > 
2m^+ + 771^^0 since the lowest threshold channel is to three pions: XiXi — ^ — > 7r~^iT^7r^. Also 
note that by generating positrons through muon decays rather than directly allows gamma ray 
constraints from final state radiation [12] to be easily evaded. 

It has been shown that a ~100 MeV dark matter particle annihilating through an a-term 
(low velocity) cross section can simultaneously yield the measured relic density and generate 
the number of positrons needed to accommodate the INTEGRAL/SPI data [10]. These are 
precisely the features of a 106-207 MeV neutralino combined with the presence of a 100 MeV-1 
GeV CP-odd Higgs. 

The main difficulty with this scenario comes from the constraints on Upsilon decays, which 
we discussed in Sec. VE. To evade the CLEO limit [43] of BR{T 7A1) < 2 x 10"^ in 
this mass region, we must require cos"^ 6a tan"^ /3 < 0.13 [see Eq. (30)]. Given these constraints, 
and considering a bino-like neutralino with a 6% higgsino admixture and m~o = 150 MeV, the 
annihilation cross section needed to avoid overproducing dark matter can only be attained for a 
fairly narrow range of rriAi ~ 2m~()±10 MeV. This scenario, although not particularly attractive 
due to this requirement, does demonstrate that it is possible to generate the INTEGRAL signal 
with neutralinos in the NMSSM. This can be confirmed or ruled out by improving the limit 
on BR{T 7^1) where the Ai is not observed or where the Ai decays to a muon pair. In 
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the latter case, the Ai may have a significant displaced vertex of a few cm, especially for small 
tan/3 and < 2m~o [48]. 

An Ai this light (300 MeV) is too light to be technically natural, however. Radiative 
corrections pull up its mass and a cancellation between different orders in perturbation theory 
is required for Ai to be this light. While we have found parameter points capable of yielding the 
INTEGRAL signal, we find that they are not stable in the sense that if any of the Higgs-sector 
parameters are adjusted by a very small amount, the Ai is pulled up in mass to (9(10 GeV). 
From our numeric analysis, rriAi as small as a few GeV is technically natural. 

X. CONCLUSIONS 

In this article, we have studied the possibility of light neutralinos (100 MeV to 20 GeV) being 
present within the Next-to-Minimal-Supersymmetric-Standard-Model (NMSSM) without con- 
flicting with constraints on the dark matter of the universe. We find that light CP-odd Higgs 
bosons with a substantial non-singlet component, which appear naturally within this frame- 
work, can provide an efficient annihilation channel for light, bino or singlino-like neutralinos. 
This channel makes it possible for very light neutralinos to generate the observed dark matter 
abundance, unlike in the case of neutralinos in the MSSM. 

Within this model, we have discussed the implications of light neutralinos for direct detection 
and find that the NMSSM can naturally provide neutralinos in the mass range (6-9 GeV) as 
required to reconcile the DAMA claim of discovery with the limits placed by CDMS and other 
experiments. We have also explored the possibility that the 511 keV emission observed from 
the galactic bulge by INTEGRAL/SPI could be generated through neutralino annihilations 
into muon pairs. This scenario appears possible for a very light (106-207 MeV) neutralino and 
for a light CP-odd Higgs boson with mass close to twice the neutralino mass, provided tan j3 is 
not large. 

This kind of scenario containing a light neutralino and/or light axion-like particles represents 
a challenge for the LHC and ILC, and is deserving of further analysis. We note that only the 
ILC will be able to study the properties of the Xi ^i-nd Ai adequately to verify that they are 
consistent with the observed dark matter density. 
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